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Abstract: Recent BaBar limits on BR(T(3S) -> ja -> jt + t~) and BR(T(3S) -> 7a 
7/Lt + /x _ ) provide increased constraints on the abb coupling of a CP-odd Higgs boson, a, 
with m a < Mx(3S). We extract these limits from the BaBar data and compare to the 
limits previously obtained using other data sets, especially the CLEO-III BR(T(1S) — > 
7 — > t + t~) limits. Comparisons are made to predictions in the context of "ideal"-Higgs 
NMSSM scenarios, in which the lightest CP-even Higgs boson, hi, can have mass below 
105 GeV (as preferred by precision electroweak data) and yet can escape old LEP limits by 
virtue of decays to a pair of the lightest CP-odd Higgs bosons, hi — > aiai, with m ai < 2niB- 
Most such scenarios with m ai < 2m T are eliminated, but the bulk of the m ai > 7.5 GeV 
scenarios, which are theoretically the most favored, survive. We also outline the impact of 
the new ALEPH LEP results in the e+e" -» Z + 4r channel. For tan/3 > 3, only NMSSM 
ideal Higgs scenarios with > 98 GeV and m ai close to 2mg satisfy the ALEPH limits. 
For tan/3 < 2, the ALEPH limits are easily satisfied for the most theoretically preferred 
NMSSM scenarios, which are those with m ai close to 2ms and m/ ll ~ 90 GeV — 100 GeV. 
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1. Introduction 

Many motivations for the existence of a light CP-odd o Higgs boson have emerged in 
a variety of contexts in recent years. Of particular interest is the m a < 2m,B region, for 
which a light Higgs, h, with SM-like WW, ZZ and fermionic couplings can have mass rrih ~ 
100 GeV while still being consistent with published LEP data by virtue of h — > 00 — )• 4r 
or 4 jet decays being dominant [Q, ||, [3|, Q (see also || |6|). Such a light Higgs provides 
perfect agreement with the rather compelling precision electroweak constraints, and for 
BR(h — > aa) > 0.75 also provides an explanation for the ~ 2.3a excess observed at LEP 
in e + e~ —> Zbb in the region ~ 100 GeV. This is sometimes referred to as the "ideal" 
Higgs scenario. More generally, superstring modeling suggests the possibility of many light 
a's. In this note, we update the analysis of J5J (see also M), quantifying the increased 
constraints on a general CP-odd a arising from recent BaBar limits on the branching ratio 



for T(35) — > 7a — > ^t + t^ decays |§ and T(3S) — > 7a — >■ decays [|10[. We also 

quantify the impact of these constraints, as well as the impact of the new ALEPH LEP 
results in the e + e~ — ?• Z + 4r final state [|ll[], on the Next-to-Minimal Super symmetric 
Model (NMSSM) ideal Higgs scenarios. 

The possibilities for discovery of an a and limits on the a are phrased in terms of the 
afi~n + , ar~r + , abb and ati couplings defined via 

^/7"^//S^ 775/a - (L1) 

(Note: when discussing a generic CP-even (CP-odd) Higgs boson, we will use the notation 
h (a). When specializing to the NMSSM context, we will use h\,h2,h^ (01,02) for the 
mass ordered Higgs states.) In this paper, we assume a Higgs model in which C „-„+ = 



- 1 - 



C aT ~ T + = C abb , as typified by a two-Higgs-doublet model (2HDM) of either type-I or type- 
II, or more generally if the lepton and down-type quark masses are generated by the same 
combination of Higgs fields. However, one should keep in mind that there are models in 
which r = (C a „-„+ = C aT - T +) / C abb S> 1 — such models include those in which the muon 
and tau masses are generated by different Higgs fields than the b mass. In a 2HDM of 
type-II and in the MSSM, C aA1 - M + = C ar - T + = C abb = tan/3 (where tan/3 = h u /hd is the 
ratio of the vacuum expectation values for the doublets giving mass to up-type quarks vs. 
down-type quarks) and C at j = cot /3. These results are modified in the NMSSM (see, e.g. 



1J and pj ). 1 In the NMSSM, both C aitl and C aibb = C aifl -^+ = C aiT - T + are multiplied 



by a factor cos6>a, where cos 6 a is defined by 

«i = cos A a M ssM + sin 6 A as, (1-2) 

where a\ is the lightest of the 2 CP-odd scalars in the model. Above, cimssm is the CP-odd 
(doublet) scalar in the MSSM sector of the NMSSM and as is the additional CP-odd singlet 
scalar of the NMSSM. In terms of cos^, = C aiT - T + = C aibb = cos 6a tan j3 and 

C ait j = cos 6a cot j3. Quite small values of cos 6a are natural when m ai is small as a result 
of being close to the U{1)r limit of the model. In the most general Higgs model, C^-^+j 
C aT - T + , C abb and C at j will be more complicated functions of the vevs of the Higgs fields and 
the structure of the Yukawa couplings. In this paper, we assume C ail -^+ = C aT - T + = C abb 
and C abl/ C att =tan 2 /3. 

For the analysis presented in this paper, we neglect the possible presence of large 
corrections at large tan/3 to C abb from SUSY loops [16, [l?], |lS[| . These are typically char- 
acterized by the quantity which is crudely of order 16 ^> ^/"^ • The correction to the 
coupling takes the form of 1/(1 + A&). Since fi can have either sign, C abb can be either 
enhanced or suppressed relative to equality with C aT - T + (the corrections to which are much 
smaller) and C a/J -^+ (the corrections to which are negligible). This same correction factor 
would apply to C aibb in the NMSSM case. 

Key ingredients in understanding current limits are the branching ratios for a — > t + t~ 
and a — > fi + decays. These branching ratios are plotted in Figs, [l] and ||. (It is important 
to note that at tree- level the a branching ratios apply equally to the a\, independent of 
COS6U, due to the absence of tree-level a,a\ — > VV couplings and similar.) Note that 
BR{a — > t + t~) and BR{a — > /i + /x~) change very little with increasing tan/3 at any given 
m a once tan/3 > 2. We note that in the region m a < 2m T , BR(a — > fi + fi~) has some 
significant structures that arise from the fact that BR(a — > gg) is substantial and varies 
rapidly in that region. The rapid variation in BR(a — > gg) occurs when m a crosses the 
internal quark loop thresholds. At higher m a , BR(a — > gg) becomes significant for m a 
near 2m},. We plot BR(a — > gg) in Fig. ||. Note that in the calculation of BR(a — > gg) 
we have chosen to keep the loop quark masses equal to the current quark masses in our 
calculations, whereas we employ thresholds of 2m^ and 2mr> for the strange quark and 
charm quark final states, respectively. Some changes in the structures present, especially 
in BR{a — > /x + /i _ ), take place if, instead, the loop quark masses are set equal to the true 
physical threshold masses. 



1 A convenient program for exploring the NMSSM Higgs sector is NMHDECAY JlJ, 
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Figure 1: BR(a — > t + t~) is plotted as a function of m a for a variety of tan/3 values. BR(a — > 
t + t~) is independent of cos 9a- 
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Figure 2: BR(a — > ) is plotted as a function of m a for a variety of tan/3 values. BR(a — > 

is independent of cos^. 



Of course, the above branching ratios are impacted by the a — > cc and a — > ss channels, 
the latter being a rather important competitor for smaller tan f3 and m a > 2mx ■ Plots of 
these branching ratios appear in Figs. || and ||, respectively. 

It is relevant to note that both BR{a — > /i + /U~) and BR{a — > r + r~) tend to decline 
slowly as m a is increased, with a significant dip in the latter for m a close to Ira^ where the b- 
loop contribution to the gga coupling is close to the point at which the internal 6's can go on- 
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Figure 3: BR(a — >■ gg) is plotted as a function of m a for a variety of tan/3 values. 
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Figure 4: BR(a — > cc) is plotted as a function of m a for a variety of tan j3 values. 



shell . This has important implications for using these channels to probe the 9 GeV < m a < 
2m b region in which many parameter choices lead to absence of light-ai finetuning in the 
NMSSM. "Light-Oi" finetuning is characterized numerically by a quantity we call G, defined 
in ||, that gives the degree of precision with which the A\ and A K soft-SUSY-breaking 
NMSSM parameters must be chosen in order that m ai < 2mg and BR(h\ — > a\a\) > 0.75 
as required to allow < 105 GeV to be consistent with published LEP constraints when 
the h\ has SM-like h\ZZ coupling. Absence of light-ai finetuning is equivalent to G <J 20. 
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Figure 5: BR(a —> ss) is plotted as a function of m a for a variety of tan/3 values. 



Typically, this condition is satisfied only when the light a\ of the NMSSM is mainly singlet. 
For example, at tan/3 = 10, 0.6 < < 1.2 (0.06 <, \ cos8 A \ < 0.12) is required if 

G < 20 is imposed as well as requiring m ai < 2ms and BR{h\ — > aia{) > 0.75, with G < 10 
achieved only for cos#a G [—0.08,-0.1], corresponding to |C b ^| G [0.8,1]. The G < 10 
range for tan/3 = 3 is broader, cos 6a G [—0.28, —0.08], while that for tan/3 = 50 is narrow, 
cos9 A G [-0.04,-0.06], yielding \C aibl \ G [0.24,0.84] and |C ai J G [2,3], respectively. 
Thus, lower tan/3 values will be harder to probe using direct limits on the a\. 

We emphasize that, given the importance of the exact a or a\ branching ratios in 
the analyses that follow, additional attention to the most precise predictions possible is 
warranted. Our a, a\ decay results employ a branching ratio program that is taken from 



HDECAY [IS]. We note that the a\ branching ratios obtained using this program are 
somewhat different than those that one obtains using the a\ decay formulae in the current 
version of NMHDECAY. In particular, the former often predicts smaller BR(a\ — > t + t~) 
than does the latter. 



2. Upsilon decay limits compared to NMSSM predictions 

Before continuing with the general analysis, it is useful to compare the limits of || and [10] 



with the predictions of the NMSSM. This comparison is done for the same two types of scans 
as in the earlier paper |2(J, except that here we focus on the 35 state rather than the 15 
state. In both scans, we hold the gaugino soft-SUSY-breaking parameters of the NMSSM 
fixed at M\ y 2${yriz) = 100,200,300 GeV and fix tan/3. In the first type of scan, called a 
"fixed-// scan", we scan over the NMSSM soft-SUSY-breaking Higgs potential parameters 
A\ and A R keeping the effective fi parameter of the model fixed at the representative value 
of n = 150 GeV (at tan/3 = 10 and 50) or \i = 152 GeV (at tan/3 = 3 for which we 
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Figure 6: BR(T(3S) - 
m ai '■ medium grey (red) 



7<zi) x BR(ai t + t _ ) for NMSSM scenarios with various ranges for 
= 2m T < m ai < 7.5 GeV; light grey (green) = 7.5 GeV < m ai < 8.8 GeV; 
and black = 8.8 GeV < m ai < 2ms GeV. The plots are for tan/3 = 3,10,50, respectively. The 
left-hand window in each plot shows results for a "fixed- //-scan" as defined in the text (and in 
Rcf. ^o|) The right-hand window shows results for F < 15 points found using a "full scan" as 
defined in the text. 
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Figure 7: We plot BR(T{3S) ->• 7 ai) x BRfa 
procedures as described in the caption of Fig. |[ 



a l non— singlet fraction (cos0 A ) 

l-i + fi~) using the same notation and scanning 



must take /i = 152 GeV in order to get physically allowable scenarios). In addition, in 
the fixed-// scans we have kept the scalar soft-SUSY-breaking masses fixed at common 
value of Msusy = 300 GeV and the A soft-SUSY-breaking parameters fixed to a common 
value of —300 GeV. In the second type of scan, termed a "full scan", we have allowed 
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/x to vary and have also allowed the soft-SUSY-breaking scalar masses and A parameters 
to vary (independently of one another). In the full scan results presented we have kept 
only scenarios with very low electroweak finetuning, as characterized by the parameter F 
(see Q for more details) being smaller than 15, where F < 15 corresponds to absence 
of electroweak finetuning. F < 15 scenarios only arise for Ss 105 GeV and are thus 
automatically "ideal" in the precision electroweak sense. As part of the fixed-/u scans and 
the full scans, we have required that the CP-even h\ escape published LEP limits by virtue 
of dominant h\ — > a\a\ — > 4r or 4 jet decays. In the forthcoming plots, the left-hand 
windows correspond to fixed-/U scan results and the right-hand windows give the results of 
a full scan for the same tan /3 value. 

Our results for the t + t~ final state are shown in Fig. ^ and those for the fi + fi~ 
final state are shown in Fig. ^. Let us focus first on the t + t~ final state. The 90% CL 
BR(T(3S) 7a) x BR(a — > t + t~) limits from BaBar range from ~ 10~ 5 at m a just 
above 2m T with a long plateau at the 3 — 7 x 10 -5 until m a passes above 10 GeV where the 
limit is of order 10~ 4 . In Fig. |(| the black points have high m ai (8.8 GeV < m ai < 2ms), 
the light grey (green) points have 7.5 GeV < m ai < 8.8 GeV and the medium grey (red) 
points have 2m T < m ai < 7.5 GeV. Let us first discuss tan/3 = 10 results, since these can 
be compared to those for Y(IS') — > 701 — > jt + t~ presented in Ref. |20f| . From comparing 
the BaBar limits summarized above with the relevant plot of Fig. ^, we see that most of 
the m ai < 7.5 GeV points are excluded, about half of the 7.5 GeV < m ai < 8.8 GeV are 
excluded, but that many fewer of the m ai > 8.8 GeV points are excluded. Still, exclusions 
of this higher m ai region are much superior to those from the CLEO-III Y(IS') data pi|| , 
which excluded none of the black points, a small fraction of the green points and about 
half of the red points. This ability to probe to higher m ai using the T(3S) is particularly 
relevant in the NMSSM context since the GUT-scale tunings of A\ and A K needed to obtain 
m ai < 2ms while at the same time having BR(h\ — > a\ai) J> 0.7, as required in the ideal 
Higgs scenario, is minimal for m ai values close to 2itib- For tan/3 = 50, one finds that 
almost all the 2m T < m ai < 8.8 GeV scenarios are excluded, but that lots of m ai > 8.8 GeV 
points survive. In contrast, for tan /3 = 3 the BaBar results only significantly constrain the 
region 2m T < m ai < 7.5 GeV. 

We now turn to the h + /j,- final state. The 90% CL BR(T(3S) -> 7a) x BR(a ->■ 
limits from BaBar are ~ 1 — 3.5 x 10 -6 for m a < 1 GeV, ~ 1 — 2 x 10~~ 6 for 1 < m a < 2m T , 
~ 1 -3 x 10~ 6 for 2m T <>m a <, 7.5 GeV, and ~ 1 - 5 x 1(T 6 for 7.5 GeV < m ai < 9.2 GeV. 
In Fig. ^ the black points have high m ai (8.8 GeV < m ai < 2m b), the light grey (green) 
points have 7.5 GeV < m ai < 8.8 GeV, the medium grey (red) points have 2m T < m ai < 
7.5 GeV and the darker grey (blue) points have m ai < 2m T . At tan/3 = 3, the final 
state data eliminates more than 4/5 of the NMSSM model points in the m ai < 2m T mass 
range, but only a small number of the NMSSM points for 2m T < m ai < 7.5 GeV and none 
of the points with 7.5 GeV < m ai . At tan/3 = 10, all m ai < 2m T NMSSM points are 
eliminated by the data as well as a small fraction of the 2m T < m ai < 7.5 GeV and 

7.5 GeV < m ai < 8.8 GeV points, but none of the 8.8 GeV < m ai points. At tan [3 = 50, all 
m ai < 2m T NMSSM points are again eliminated, perhaps half of the 2m T < m ai < 7.5 GeV 
points are eliminated, a still significant fraction of the 7.5 GeV < m ai < 8.8 GeV points 
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are eliminated, and even a significant number of the 8.8 GeV < m ai points are eliminated. 

To summarize, only the /x + pT channel provides constraints for m ai < 2m T and almost 
all the ideal-Higgs-like NMSSM scenarios with tan/3 > 3 are eliminated. For 2m T < m ai , 
the t + t~ channel provides the most eliminations for all tan/3. Certainly, the BaBar T(35) 
results are a big stride relative to the CLEO-III Y(IS') results, especially at m ai < 2m T and 
at high m ai . Of course, it is important to note that the NMSSM scenarios most favored in 
order to minimize light-ai finetuning have m ai very near 2m# and thus cannot be limited 
by Upsilon decays. 



3. General limits on the abb coupling 

Our ultimate goal is to use the T35 limits in combination with other available limits to 
extract limits on \C bb \. The older experiments that provide the most useful constraints 
are as follows. Prior to the recent BaBar data, for 2m T < m a < 9.2 GeV the recent CLEO- 



III [21] limits on T(1S) — > 7a — > jt + t were the strongest. For 9.2 GeV < m a < Mwis), 



mixing of the a with various 1% and xo bound states becomes crucial [22]. Ref. [21] gives 



results for C™| x in this m a range without taking this mixing into account but notes that 
their limits cannot be relied upon for m a > 9.2 GeV. Whether additional limits can be 
extracted from lepton non-universality studies in the 9.2 < m a < Mxns) region is being 
studied Q. OPAL limits @ (which assume BR(a -s- t+t~) = 1) on e + e~ -> bbr + T~ 



become numerically relevant for roughly 9 GeV < m a < 2m b- Ref. [24] converts these 



limits to limits on the abb coupling using the modeling of [22]. These are the only LEP 
limits in the Myr^g) < m a < 2rriB range and continue to be relevant up to 12 GeV. 
Above m a = 2rriB these abb coupling limits become quite weak due to the % — a mixing 
uncertainties and the decrease of BR(a —¥ t + t~). For m a > 12 GeV, limits on the abb 



coupling can be extracted from e + e~ — >■ bba — > bbbb [25]. One should also keep in mind 
that values of |C ab ^| above 50 raise issues of non-perturbativity of the abb coupling and 
are likely to be in conflict with Tevatron limits on bba production |2q] . The limits, C^| X ! 
on C h t coming from all data, including the recent BaBar results, are plotted in Fig. H 



for various Ry t = \/C ab -^/C at j values. (In a 2HDM model type-II context, R b u = tan/3.) 
Note the rapid deterioration as m a — > My (35). The variation with R b n arises because 
BR(a — >■ t + t~) varies with R b u as shown in Fig. |l]. Basically, for tan/3 > 1 the BaBar 
results provide the most stringent limits. For tan/3 = 0.5 the a decays to a complicated 
mix of channels and the old CUSB-II limits (which were independent of the exact a final 
state) are strongest for m a < 8 GeV. 

In Fig. H|, we have also plotted limits extracted p^] from Tevatron data using a reinter- 



pretation of a CDF analysis performed over the range 6.3 GeV < M^+^- < 9 GeV [g8|, |29]. 
This analysis placed limits on the ratio R = /-J^ £ \ppfc~^ ^J--, , where e was a narrow 
resonance produced in the same manner as the Tig. Fluctuations of R above a smooth fit 
to the overall spectrum were searched for and 90% CL limits were placed on R. It is rel- 
atively straightforward to apply this analysis to place limits on R = g (Yi°)Bfl(Tis->-/j+^-) " 
The 90% CL limits on it! corresponding to the available L = 630 pb _1 data set are then 
easily converted to limits on \C abb \- These limits as a function of m a are those plotted as 
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Figure 8: Upper limit, C™| x , on |C ab j| as a function of m a for a variety of tan/3 values coming 
directly from experimental data. The highest (red) curve is for tan/3 = 0.5, the other curves, in 
order of decreasing are for tan/3 = 1, tan/3 = 2 and tan/3 > 3. 

° abb — 

the solid histogram. A simple statistical extrapolation of these limits to L = 10 fb -1 (an 
integrated luminosity that will soon be available) is shown as the dashed histogram. These 
limits hold for tan f3 > 2. We see that these limits improve rapidly as m a increases. While 
the L = 630 pb _1 limits are not quite competitive with the limits from BaBar data at 
m a ~ 9 GeV, we observe that the L = 10 fb _1 limits will actually be slightly better if the 
extrapolation holds. 

While T(n/S)-based limits are kinematically limited and become weak for m a > 9.6 GeV, 
there is no such kinematic limitation for limits based on hadronic collider data. In fact, CDF 
measured the M /J + /J - spectrum above 9 GeV, but did not perform the easily reinterpreted 
R analysis in the region M At + At - > 9 GeV. In |27|, we estimated the 90% CL limits from 
the L = 630 pb _1 measurements in the M„+„- > 9 GeV region (out to M^+ M - = 12 GeV) 
and found that, in the range 9.6 GeV m a <^ 2mg, implied limits on |C a ^| were of or- 
der \C ab i\ < 1.6 — 1.8 for m a outside the Y 2 s and T35 peaks. At both peaks we found 
\C abl \ < 2. For L = 10 fb~\ these limits should come down to \C ab j ) \ < 1, and begin to 
constrain the most preferred NMSSM parameter regions, especially for large tan/3. 

4. Implications of general abb limits for NMSSM scenarios 

In the NMSSM, we note that it is always possible to choose cos 6a so that the limits on 
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Figure 9: cos6*™ ax in the NMSSM (where C ab j = cos 9a tan /3) as a function of m a . The different 
curves correspond to tan/3 = 1 (upper curve), 1.7, 3, 10, 32 and 50 (lowest curve). CDF/Tevatron 
constraints do not affect this plot. 



C ai bb as a function °f tan/3 are satisfied. The maximum allowed value of | cos 9a\, cos #™ ax , 
as a function of m a = m ai for various tan /3 values is plotted in Fig. [| Constraints are 
strongest for m a < 9 GeV for which Upsilon limits are relevant, and deteriorate rapidly 
above that. As seen in Fig. |||, currently the limits from the Tevatron/CDF data are not as 
strong as those from the BaBar T35 data and do not affect this plot. 

As an aside regarding the general 2HDM(II) model, we note that any point for which 
cos #™ ax is smaller than 1 corresponds to an m a and tan f3 choice that is not consistent with 
the experimental limits. Disallowed regions emerge in the range m a <> 2m T for tan/3 = 1, 
rising quickly to m a <, 9 GeV for tan/3 = 1.7 and m a <> 10 GeV for tan/3 > 3. These 
excluded regions apply to any light doublet CP-odd Higgs boson, including the beyond the 
MSSM scenarios of [H], 32] which are consistent with other experimental constraints 
for tan j3 ^ 2.5. 

We can illustrate the effects of the limits plotted in Fig. |9| on preferred NMSSM sce- 
narios. Relevant plots appear below. The first set of plots, Figs. [l0|, [II] and 12, for 
tan/3 = 3, 10, and 50, respectively, show results for "fixed-/i scans" (see earlier definition). 
In each figure, the left-hand plot gives the light-ai finetuning measure G as a function 
of cos 9a before imposing the cos # max constraint while the right-hand plot gives G as a 
function of cos 9 a after imposing cos# max . The point notation is according to m ai : blue 
for m 0l < 2m T , red for 2m T < m ai < 7.5 GeV, green for 7.5 GeV < m ai < 8.8 GeV and 
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Figure 10: Light-ai fmetuning measure G before and after imposing limits |cos#a| < cos 9™^. 
These plots are those obtained for "fixed-^ scans" with /i = 152 GeV and setting tan/3 = 3. 
Note that many points with low m ai and large | cos6U| are eliminated by the | cos 8 a \ < cosO™ 3 * 
requirement, including almost all the m ai < 2m T (blue) points and a good fraction of the 2m r < 
m ai < 7.5 GeV (red) points. 



tan/3=10, fj,= 150 GeV, M 123 = 100,200,300 GeV 




tan/S=10, /i=150 GeV, M 123 = 100,200,300 GeV 




Figure 11: As in Fig. [l0| but for (i = 150 GeV and tan/3 = 10. Note that many points with 
low m ai and large | cos 6a \ are eliminated, including almost all the m ai < 2m T (blue) points and 
2m T < m ai < 7.5 GeV (red) points. 



black for 8.8 GeV < m ai < 2ms ■ We see that the bulk of points with m ai < 7.5 GeV are 
eliminated by the costf™ 331 " limit and that the points with m ai > 7.5 GeV at large | cos 9a\ 
are also eliminated. 



The second set of plots below, Figs. O, 14 and Qq, show results for "full scans" , as 



defined previously, for tan /3 = 3, 10, and 50, respectively. Only points with electroweak 
fmetuning measure F below 15 are plotted. As in the previous set of plots, the left-hand 
plot in each figure shows the points allowed without the cos #™ ax constraint and the right- 
hand plot displays the points remaining after imposing cos #™ ax . The limited statistics for 
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Figure 12: As in Fig. [j^, but for /x = 150 GeV and tan/3 = 50. Note that the only surviving 
points are those with m ai > 8.8 GeV (black points) at small | cos 9a\ ■ 
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Figure 13: Light-ai finetuning measure G before and after imposing | cos9a\ < cos0™ ax . These are 
the results obtained using a "full scan" at tan fj = 3. Only solutions with electroweak finetuning 
measure F < 15 are retained. Note that a good fraction of the m ai < 2m T (blue) points and 
2m T < m 0l < 7.5 GeV (red) points are eliminated by the cos#™ ax cut. 

the parameter scans that search for points with low F are apparent, but the trends are 
clearly the same as in the fixed \i scans presented previously. 

From a theoretical perspective, an interesting pattern emerges: the cos #™ ax constraint 
eliminates those points for which the light-ai finetuning measure G is never small and 
zeroes in on those cos 6a values for which small G is quite likely. 

5. Effective £ 2 in the h — > 4r channel for vector-boson fusion at the LHC 
and LEP Zh channel constraints 

Discovery of a Higgs using vector boson fusion at the LHC or at LEP with 2m T < m ai < 
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Figure 14: As in Fig. [T|, but for tan/3 = 10. Note that many points with lower m ai and large 
| cos6U| are eliminated by the | cos 6 a | < cos#™ ax cut. 
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Figure 15: As in Fig. but for tan/3 = 50. Note that no F < 15 points found in our scans 
survive the | cos 6 a \ < cos#™ ax limits. 

2m b (which is the only kind of point that survives with G < 20) is essentially determined 
by 

? = fife) BR(h -> aa)[BR{a -► t+t~)] 2 . (5.1) 
\9vv / 

We consider expectations for £ 2 in the NMSSM ideal Higgs scenarios with the cos#™ ax 
constraint imposed in addition to the usual constraints contained within NMHDECAY. 

In Fig. [16| we take tan (5 = 3 and plot £ 2 for h = h\ and a = a% as a function of m 01 
and as a function of for points coming from the fixed fi scans after imposing G < 20 
and requiring | cos9a\ < cos#™ ax (m a ). We observe that £ 2 as small as ~ 0.42 is possible at 
high m ai , which points tend to have m/ ll € [90, 100] GeV. As seen in Fig. [I?], these same 
remarks apply also to the F < 15 points obtained in our finetuning scans when G < 20 and 
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Figure 16: £ 2 for h = hi as a function of m ai and for points with G < 20 and | cos Qa\ < 
cos#™ ax (m a ). These plots are those obtained using the "fixed- fi" scanning procedure for tan/? = 3. 
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Figure 17: £ 2 for h = hi as a function of m ai and for points with F < 15, G < 20 and 
| cos 9a\ < cos0™ ax (m a ). These plots are those obtained using the described scanning procedure 
for tan/3 = 3. 

| cos 9a\ < cos#™ ax (m ai ) are imposed. These same remarks also apply to the tan/3 = 10 
plots of Figs. |l8| and [l9| as well as to the tan /3 = 50 fixed- //-scan plot of Fig. |2^. (Note 
that no F < 15, G < 20 points survived our limited statistics electroweak finetuning scan 
in the tan /3 = 50 case and so there is no corresponding figure.) 

In addition, we have also considered £ 2 expectations in scenarios with rather low tan/3. 
These were detailed in |3^]. There, we performed fixed-// scans as defined earlier, with the 
difference that at tan/3 = 1.7 and tan/3 = 1.2 we used different values for M$usy and 
A parameters, which values are indicated on the figures. At tan /3 = 2 we employed 
Msusy = —A- = 300 GeV as for the fixed-// scans for tan/3 = 3, 10, 50. 

The main distinguishing characteristic of the low tan /3 scenarios is that both h\ and 
hi can be light with masses not far from 100 GeV, although there are certainly choices for 
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Figure 18: £ 2 for h = hi as a function of m ai and m/^ for points with G < 20 and | cos 6a\ < 
cos 6*™ ax (m a ). These plots are those obtained using the "fixed- /i" scanning procedure for tan (3 = 10. 
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Figure 19: £ 2 for h = hi as a function of m ai and for points with F < 15, G < 20 and 
| cos 0a | < cos^™(m a ). These plots are those obtained using a "full scan" for tan/3 = 10. 



the NMSSM parameters for which only h\ is light while h,2 is much heavier. When /12 is 
light, the charged Higgs can also have mass close to 100 GeV. 2 Here, our interest is in 
the predictions for £ 2 . 

Results for £ 2 at tan /3 = 2 are rather similar to those found for higher tan /3, as shown in 



Fig. 21|. In this figure, the blue +'s are all points that satisfy the NMHDECAY constraints 
- unlike the previous figures, color coding is not employed to distinguish different m ai 
values. Results for £2 are n °t shown; even when m/ 12 is close to 100 GeV, £ 2 is quite small. 
This tan/3 = 2 case is similar to the tan/3 = 3, 10, 50 cases also in that it is almost always 



2 Note that a light H can cause the NMSSM prediction for BR(b — > sj) to substantially exceed the 
experimental value, which is only slightly above the SM value. Thus, contributions from other SUSY 
diagrams must enter to cancel the H diagrams. In models with low finetuning, SUSY is light and such 
cancellation is generically entirely possible. 
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Figure 20: £ 2 for h = hi as a function of m ai and m/^ for points with G < 20 and | cos Qa\ < 
cos 6*™ ax (m a ). These plots are those obtained using the "fixed-/i" scanning procedure for tan j3 = 50. 
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Figure 21: £f as a function of m ai and for points with G < 20 and | cos 0a | < cos6™ ax (m ai ) 
and tan/3 = 2. These plots are those obtained using a "fixed-//' scanning procedure with the [a, 
Msusy and A parameters indicated on the figure. We have not indicated different m ai mass ranges 
using different colors in these figures. 



the case that VV couples primarily to the h\ so that when m/ ll < 105 GeV we have the 
"ideal" Higgs explanation of the precision electroweak data. 

For tan/3 <J, 1.7, there are some interesting new subtleties compared to tan/3 J> 2. 
Plots of £f of the h\ and £| °f ^2 appear in Figs, 
we follow the notation established in Ref. [p3| ]. In detail, the blue +'s are all points that 
satisfy the NMHDECAY constraints. The red crosses single out those points for which 

< 65 GeV. Yellow squares indicate points for which BR{h\ — > a\a\) < 0.7. In p3| . 
there were also points indicated by green diamonds for which in addition the light CP-odd 
Higgs is primarily doublet-like, cos 2 9a > 0.5. However, these are absent from the present 
plots, not because of the improved cos #™ ax limits from the recent BaBar data, but rather 
because of the G < 20 requirement which very strongly disfavors large | cos 8a\ at all m ai , 



22 and E3, respectively. In these plots, 
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Figure 22: £f as a function of m ai and mi ll for tan/3 = 1.7 points obtained from a nxed-^t scan 
after requiring G < 20, m e // < 100 GeV and |cos#a| < cos 6*™ ax (m ai ). The point notation is 
explained in the text. 



including m ai above Myr 3 gy Of course, the BaBar data eliminates many points with 
m ai < M T (35) having cos 2 9a < 0.5, certainly more than in the analysis of 

Let us now discuss the tan/3 = 1.7 case in more detail. We first wish to discuss 
the extent to which the points that survive the NMHDECAY scans can be "ideal" in the 
precision electroweak sense. Defining 

CVi = gvvhi/gvvhsM ' CV2 = gwhj gvvKsu > ( 5 - 2 ) 

then, noting that it is a good approximation to neglect any /13 coupling to VV, one has 
the sum rule 

CVf + CVi - 1 • (5-3) 

In this notation, the effective precision electroweak mass, m e ff, is given to very good 
approximation by 

CV, 2 CK 2 -s 

m e// = m ftl 1 m^ 2 2 . (5.4) 

In order to guarantee that all accepted points are ideal, we require as part of our tan (3 = 1.7 
scan that m e ff < 100 GeV. 3 Now, let us describe the associated plots. First, very low 
values of are possible (see the red crosses). These red cross points are such that £ 2 
and are comparable and both below 0.2. Second, very few of the yellow square points 
(defined by BR(h\ — > a\ai) < 0.7) survive the ideal requirement. But, those that do have 
quite small £ 2 and £|. The run-of-the-mill blue + points have somewhat larger £ 2 < 0.4 
and somewhat smaller £ 2 < 0.2. Overall, the 4r final state in h± and ti2 decays typically 
has significantly smaller cross section for tan/3 = 1.7 as compared to tan/3 > 2. 

The lowest value of tan/3 consistent with maintaining perturbativity up to the GUT 
scale is tan /3 = 1.2. £ 2 and £| plots for this case appear in Figs. |24] and 25, respectively. In 
this case, the effective £ 2 values are mostly quite small. Relative to the tan/3 = 1.7 plots, 
the main thing that has changed is that BR{a\ — > t + t _ ) has declined substantially. The 
majority of the a\ decays are into gg and cc, i.e. final states that are harder to constrain. 



3 This was not imposed in the plots of |33| . 
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Figure 23: £f as a function of m 0l and mh 2 for tan/3 = 1.7 points obtained from a fixed-^t scan 
after requiring G < 20, m e // < 100 GeV and | cosff^l < cos 0™ ax (m Ql ). 
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Figure 24: £f as a function of m ai and for tan/? = 1.2 points obtained from a fixed-^ scan 
after requiring G < 20, m e ff < 100 GeV and | cos 0a\ < cos 6>™ ax (m Ql ). 

Of course, the knowledgeable reader will recognize that all the £ 2 plots presented are 
aimed at comparing these NMSSM models to the new ALEPH analysis of the 4r final 
state 0. According to the ALEPH analysis, to have m hl < 100 GeV, £f < 0.52 (0.42) is 
required if m ai ~ 10 GeV (4 GeV). These limits rise rapidly with increasing m/ ll — for 
m/ ll = 105 GeV (the rough upper limit on mh 1 such that electroweak finetuning remains 
quite small and precision electroweak constraints are fully satisfied) the ALEPH analysis 
requires £ 2 <^ 0.85 (<^ 0.7) at m ai ~ 10 GeV (4 GeV). These limits are such that the easily 
viable NMSSM scenarios are ones: i) with m ai below but fairly close to 2ms, which is, 
in any case, strongly preferred by minimizing the light-ai finetuning measure G; and/or 
ii) with tan/3 relatively small (< 2). 4 These are also the scenarios for which Upsilon 
constraints are either weak or absent. In particular, we note the following: a) all tan/3 < 2 
cases provide < 100 GeV scenarios that escape the ALEPH limits; b) there are a few 



4 A similar conclusion applies to models beyond the MSSM with a light doublet CP-odd Higgs boson |3(i| , 
gj , [32j , Since these scenarios are consistent with other experimental limits only for tan/3 < 2.5, the new 
preliminary Aleph limits only constrain the upper range of the allowed region of tan /?. 
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Figure 25: £f as a function of m Ql and m/, s for tan/3 = 1.2 points with G < 20, m e // < 100 GeV 
and |cos#a| < cosf?™ ax (m ai ). 

G < 20, tan /3 = 3 scenarios with m/ ll as large as 98 GeV and 99 GeV and with £ 2 essentially 
equal to the ALEPH limits of t; 2 < 0.42 and £ 2 < 0.45 applicable at these respective 
values; c) tan/3 = 10 ideal scenarios easily allow for ra\ tl ~ 100 — 105 GeV (because the 
tree- level Higgs mass is larger at tan/3 = 10 than at tan/3 = 3) and at m ai ^ 2niB 
many J> 100 GeV points have £ 2 < 0.5 in the fixed-// scan and a few of the full-scan 
points have £ 2 < 0.6 for ~ 105 GeV, both of which are below the m ai = 10 GeV 
ALEPH upper limits on £ 2 of 0.52 at m hl ~ 100 GeV and 0.85 at m hl = 105 GeV; d) at 
tan/3 = 50 there are some G < 20 points with ~ 100 GeV and m ai < 2mB having £ 2 
below the 0.52 ALEPH limit. Finally, we note that for the entire range of Higgs masses 
studied the ALEPH limits were actually ~ 2a stronger than expected. Thus, it is not 
completely unreasonable to consider the possibility that the weaker expected limits should 
be employed. These weaker limits for example allow £ 2 as large as 0.52 at ~ 95 GeV 
and 0.9 for ~ 100 GeV. These weaker limits allow ample room for the majority of the 
m ai <> 2m b ideal Higgs scenarios. 

6. Conclusions 

In this paper, we have updated the constraints on the NMSSM ideal Higgs scenarios in 
which h\ (and for low tan/3, also possibly /12) has mass <i 105 GeV and decays largely 
(but not entirely) via hi — > aiai — > t + t~t + t~ . Such low mass(es) for the Higgs boson(s) 
with large VV coupling are strongly preferred by precision electroweak data and are also 
strongly preferred in order to minimize electroweak finetuning. Indeed, all the NMSSM 
points plotted in this paper have effective precision electroweak mass below ~ 105 GeV. 
The new data that constrains such scenarios derives from T35 — > ^yfi + fi~ and 7t + t~ decay 
data from BaBar and ALEPH studies of the e + e~ — > ZAt final state. The latter was 
employed by ALEPH to place limits as a function of and m ai on the quantity £ 2 = 
— >■ ai<2i) [BR(ai — > r + r~)] 2 . Although these new constraints are significant, 
there is still ample room for the ideal Higgs scenarios, especially if tan (5 is small and 
m ai <i 2m,B (the latter region being that for which the "light-ai" finetuning measure is 
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minimal and also BR(a\ — > t + t~) is somewhat suppressed). For tan/3 > 3, it is only 
the m ai < 2mB points that can escape the ALEPH £ 2 limits. The case of tan /3 = 3 is 
the most marginal with only a few NMSSM points with m/ il < 99 GeV (the rough upper 
limit on at tan/3 = 3) having £ 2 essentially equal to the ALEPH limit at a given 
rrifa. For tan/3 = 10, one finds scenarios with ~ 100 — 105 GeV and £ 2 ~ 0.43 
when m ai < 2m#, which £ 2 is well below the ALEPH limit of ~ 0.52 — 0.85 for such 
and m ai . At tan/3 = 50, although our scanning statistics were limited, we found 
points with m/ ll ~ 100 GeV and m ai < 2rriB having ^ 2 below the 0.52 ALEPH limit. 
(We note that the ALEPH limits are significantly stronger than the ALEPH collaboration 
was expecting. If one were to use expected limits instead then the tan /3 > 3 scenarios 
would be much less constrained.) For tan/3 <> 2, the ideal-Higgs NMSSM scenarios are 
not particularly constrained by the ALEPH limits. In particular, for tan/3 = 2,1.7,1.2 
one finds m^, < 100 GeV scenarios with £ 2 < 0.32,0.23,0.15, respectively. The lower £ 2 
values arise because these lower tan/3 values have increasingly reduced BR{a\ — > t + t~), 
which, in turn, is due to increasingly larger values of BR(ai — > gg + cc). Such £ 2 values 
are completely consistent with the ALEPH limits. 

The Tevatron and LHC discovery prospects for the Higgs bosons in the low-tan /3 
scenarios have yet to be fully analyzed. Searches for the h± and the a\ using the a\ — > t + t~ 
and a\ — > fJ- + fJ.~ decay modes will certainly become more difficult as these branching 
ratios decline with decreasing tan/3. Such search modes include: direct (vs. coming from 
hi —7- aidi) detection of the a\ at the Tevatron and LHC in the gg — > a\ — » /i + /i~ 
channel J27J; searches for gg — > h\ — > a\a\ — > t + t~t + t~, t + t~ fi + and/or [j& fj& [f~ 



at the Tevatron [34| and LHC [pq| ; and LHC detection of pp pph\ with h\ — > a\a\ — > 
t + t~t + t~ |^6| . Backgrounds in the increasingly important channels with a\ — > gg + cc 
will undoubtedly be much larger and will make discovery employing these latter a\ decay 
modes quite difficult. 

As part of the NMSSM study, we first obtained updated limits on the abb coupling (as- 
suming C ab ^ = C ar - T + = C a/J - At +) that are applicable in a wide variety of model contexts. 
The main improvements in these general limits result from recent BaBar data. 

Finally, one should not forget that the NMSSM is only the simplest model of a general 
category of SUSY models having one or more singlet scalar superfields in addition to 
the usual two-doublet scalar superfields. Such models are generically very attractive in 
that they allow for an NMSSM-like solution to the \i problem, while maintaining coupling 
constant unification and RGE electroweak symmetry breaking as in the MSSM. In addition, 
models with more than one extra singlet scalar superfield will allow one or more light Higgs 
bosons with SM-like couplings to V V (a scenario having excellent agreement with precision 
electroweak constraints and minimal electroweak finetuning) that can escape Upsilon and 
LEP limits more easily than the NMSSM by virtue of multiple decays channels of the 
Higgs-s> a k aj, . . .type. 
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